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Abstract. We consider Fabry-Perot cavity resonance in periodic stacks of 
anisotropic layers with misaligned in-plane anisotropy at the frequency close to 
a photonic band edge. We show that in-plane dielectric anisotropy can result 
in a dramatic increase in field intensity and group delay associated with the 
transmission resonance. The field enhancement appears to be proportional to 
forth degree of the number N of layers in the stack. By contrast, in common 
periodic stacks of isotropic layers, those effects are much weaker and proportional 
to TV 2 . Thus, the anisotropy allows to drastically reduce the size of the resonance 
cavity with similar performance. The key characteristic of the periodic arrays 
with the gigantic transmission resonance is that the dispersion curve ui (k) at 
the photonic band edge has the degenerate form Au ~ (Afe) 4 , rather than the 
regular form Auj ~ (AA;) 2 . This can be realized in specially arranged stacks of 
misaligned anisotropic layers. The degenerate band edge cavity resonance with 
similar outstanding properties can also be realized in a waveguide environment, 
as well as in a linear array of coupled multimode resonators, provided that certain 
symmetry conditions are in place. 
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1. Introduction 

The subject of this paper is the Fabry-Perot cavity resonance in periodic layered 
structures. This phenomenon, known for many decades, is also referred to as the 
transmission band edge resonance, because it occurs in the vicinity of photonic band 
edge frequencies in finite photonic crystals and is accompanied by sharp transmission 
peaks. This basic effect can occur in a finite periodic array of any two alternating 
materials with different refractive indices. It can also be realized in a waveguide 
environment, or in a finite array of coupled resonators. The only essential requirements 
are: (a) low absorption, (b) the appropriate number N of unit cells in the periodic 
array, and (c) the presence of a frequency gap (a stop-band) in the frequency spectrum 
of the periodic array. The latter is a universal property of almost any lossless 
spatially periodic structure. This phenomenon has found numerous and diverse 
practical applications in optics. Below we briefly outline some basic features of the 
classical band edge resonance in periodic stacks of isotropic layers. We only highlight 
those points, which are necessary for the following comparative analysis of stacks 
involving anisotropic layers. More detailed information on the relevant aspects of 
electrodynamics of layered dielectric media can be found in an extensive literature on 
the subject (see, for example, ^JElElEIElElEIBUnij an d references therein) . 

Consider a simplest periodic array of alternating dielectric layers A and B, as 
shown in Fig. ^a). The layers are made of transparent isotropic materials with 
different refractive indices ua and uq. Such periodic stacks are also referred to as 
a ID photonic crystal. The electromagnetic eigenmodes of the periodic structure in 
Fig. [Ja) are Bloch waves with a typical wave number/frequency diagram shown in 

Fig.ntb). 

Consider now a finite periodic stack composed of N unit cells L in Fig. Qfa). 
Such a stack is commonly referred to as a Fabry-Perot cavity. If the number N of the 
double layers L is significant, the stack periodicity causes coherent interference of light 
scattered by the layers interfaces. In practice, periodic stacks having as few as several 
periods L, can display almost total reflectivity at the band gap frequencies, provided 
that the refractive indices ha and ng of the adjacent layers differ significantly. 

A typical frequency dependence of the finite stack transmittance is presented 
in Fig. [21 The frequency range shown includes the vicinity of the photonic band 
edge (BE) g in Fig. ^b). The sharp transmission peaks below the photonic band 
edge frequency oj g correspond to transmission band edge resonances, also known as 
Fabry-Perot cavity resonances. At each resonance, the electromagnetic field inside the 
periodic stack is close to a standing wave composed of a forward and a backward Bloch 
eigenmodes with large and nearly equal amplitudes. The slab boundaries coincide with 
the standing wave nodes, where the forward and backward Bloch components interfere 
destructively, as illustrated in Figs. 01 and 01 The latter circumstance determines 
the wave numbers of the forward and backward Bloch components at the resonance 
frequencies 

k s~k g ±—s, a = 1,2,..., (1) 

where s is the order number of the resonant peak in Fig. [21 and k g is the wave number 
corresponding to the photonic band edge. In our case, k g = ir/L. The resonance 
frequencies themselves can be expressed in terms of the dispersion relation u> (fc) of 
the respective frequency band. Indeed, just below the photonic band edge g in Fig. 
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Figure 1. a) Periodic stack composed of two alternating layers A and B, each of 
which is made of isotropic transparent material. L is the unit cell of the periodic 
structure. The physical characteristics of the stack are specified in the Appendix, 
b) The respective k — u) diagram. The wavenumber k and the frequency u) are 
expressed in units of l/L and c/L. The point g at the Brillouin zone boundary 
designates the edge of the lowest frequency band. 



^b), the dispersion relation w (k) can be approximated by the quadratic parabola 



-f(k-k g ) , where J' g = l^ — 



< 0. 



(2) 



This relation together with Q yield the frequencies of the resonant transmission peaks 
in Fig. H 
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where tu g — to (k g ) is the band edge frequency. For example, the transmission peak 1 
closest to the photonic band edge is located at 
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The dependence (j2Jl is illustrated in Fig. |21 

Electromagnetic field distribution inside the stack at the frequency of band edge 
resonance is shown in Figs. 01 and 01 for the first two transmission resonances, 
respectively. For a given amplitude \Pj of the incident wave, the maximal field intensity 
|\t(z)| 2 inside the slab depends on the number N of the double layers in the stack 
and on the order number s of the resonance peak in Fig. |21 



■> i'N 

max|* 0)| oc |*/| ( — 



(5) 



The maximal field intensity JSJ is proportional to squared thickness of the slab and, 
for a large N, is greatly enhanced compared to that of the incident light. By contrast, 



Gigantic transmission band edge resonance in periodic stacks of anisotropic layers 4 



a) N=16 b) N = 32 




1.6 1.65 17 co 16 165 17 CO 

Frequency 9 ^9 

Frequency 



Figure 2. Typical transmission spectrum of finite periodic stacks composed of 
different number N of unit cells L shown in Fig. 1(a). The sharp transmission 
peaks below the band edge frequency ui g are associated with photonic band edge 
resonance. 

at the slab boundaries at z = and z = D = NL, the field amplitude '5 (z) always 
remains comparable to fyj to satisfy the electromagnetic boundary conditions i|47|) . 
At frequencies outside the resonance transparency peaks, the field amplitude inside 
the stack drops sharply. 

The exact definition of the physical values plotted in Figs. Oland01 as well as the 
numerical parameters of the periodic stacks used to generates these plots, are given 
in the Appendix. 

The expressions Q through ((J5J are valid if N ^> 1 and only apply to 
the transmission resonances close enough to the photonic band edge. In further 
consideration, we will focus on the most powerful first resonance s = 1, closest to 
the photonic band edge. 

Above, we outlined some basic features of the transmission band edge resonance 
in finite stacks of isotropic layers. The question we would like to address in this paper 
is whether the presence of anisotropic layers in a finite periodic stack can qualitatively 
change the nature of the Fabry-Perot cavity resonance. We will show that, indeed, 
in periodic stacks involving anisotropic layers, the transmission resonance can be 
significantly stronger, compared to what is achievable with common periodic stacks of 
isotropic layers. For instance, in the periodic stack shown in Fig. [5] the field intensity 
associated with the transmission band edge resonances can be proportional to A 4 , 
rather than N 2 . The latter implies that a stack of N anisotropic layers can perform 
as well as a common stack of N 2 isotropic layers. And this is a huge difference! The 
physical reason for this is that periodic stacks of anisotropic layers can support the 
kind of k — lu diagrams that are impossible in stacks of isotropic layers. Specifically, a 
dispersion curve u (k) of the stack in Fig. [5] can develop a degenerate band edge, as 
shown in Fig. EJb). Just below the degenerate band edge d, the dispersion curve can 
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Figure 3. Smoothed intensity distribution (59) inside periodic stack at the 
frequency u>\ (TV) in (4) of the first transmission resonance. The amplitude of 
the incident wave is unity. The distance z from the left boundary is expressed in 
units of L. At the stack boundaries at z = and z = NL, the field intensity is of 
the order of unity, while inside the slab it is enhanced by factor TV 2 . 
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Figure 4. The same as in Fig. 3, but at the frequency L02 (N) of the second 
transmission resonance in Figs. 2(a) and 2(b), respectively. 
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which implies a huge density of modes. And this is what makes all the difference 
compared to the case (0) of a regular band edge. 



Figure 5. Periodic stack capable of supporting k — ui diagram with degenerate 
band edge. A unit cell L includes three layers: two birefringent layers Ai and 
A2 with different orientations ip 1 and ip 2 of the respective anisotropy axes in the 
X — Y plane, and one isotropic B layer. The misalignment angle ip = tpj — ip 2 
between adjacent A layers must be different from and 7r/2. 

For convenience, the domain of definition of the Bloch wave number k in Fig. EJis 
chosen between and 2w/L. Since the Bloch wave number is defined up to a multiple 
of 2ir/L, the representation in Fig. EJis equivalent to that in Fig. ^b). Note that the 
points a, <?, and d in Fig. and|SJb) lie at the Brillouin zone boundary at k = tt/L. 

In Fig. \7\ we present the transmission dispersion of the finite periodic stacks in 
Fig. [S] composed of 32 unit cells L and having the k — u> diagram in Fig. HJb). The 
frequency range shown includes the degenerate band edge (DBE) at ui — ujd- The field 
intensity distribution at the frequency of the first transmission resonance 1 is shown 
in Fig. IBIb). One can see that for a given N, the resonant field intensity in Fig. [S]is 
significantly larger compared to that of the vicinity of a regular band edge, shown in 
Fig. |3 Specifically, in the case of degenerate band edge 



max ^ ( r)j 2 oc |^z| 2 ( ^ 



(7) 



compared to the estimation 101 related to a regular band edge. The transmission 
bandwidth in the case of DBE appears to be much smaller. 

In practice, the field amplitude associated with the transmission resonance is 
limited not only by the number of layers in the stack, but also by such factors as 
absorption, nonlinearity, imperfections of the periodic array, stack dimensions in the 
X — Y plane, incident radiation bandwidth, etc. All else being equal, the stack with 
degenerate band edge can have much fewer layers and, therefore, can be much thinner 
compared to a regular stack of isotropic layers with similar performance. Much smaller 
dimensions can be very attractive for a variety of practical applications. Similar 
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Figure 6. The first band of the k — ui diagram of the periodic stack in Fig. 5 for 
four different values of the B - layer thickness, Db- (a) Db = 0.75891 X L. (b) 
Db = 0.45891 X L\ in this case the upper dispersion curve develops degenerate 
band edge d. (c) Db = 0.35891 X L. (d) Db = 0; in this case the two intersecting 
dispersion curves correspond to the Bloch waves with different symmetries, and 
the respective modes are decoupled. The physical parameters of the periodic 
structure are specified in the Appendix. 
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Figure 7. The transmission spectrum of periodic stack with the k — u) diagram 
in Fig. 6(b). The stack is composed of N = 32 three-layered unit cells L in Fig. 
5. The sharp transmission peaks below the degenerate band edge frequency u>d 
are associated with Fabry-Perot cavity resonances. 
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Figure 8. Smoothed intensity distribution inside periodic stack in Fig. 5 at the 
frequency of the first Fabry-Perot resonance, closest to degenerate band edge d. 
The amplitude of the incident wave is unity. At the stack boundaries at z = 
and z = NL, the field intensity is of the order of unity, while inside the slab 
electromagnetic energy density is enhanced by factor TV 4 . 



effect associated with degenerate band edge (© of the k — u> diagram can also be 
achieved in the waveguide environment, as well as in finite periodic arrays of coupled 
multiple-mode resonators. On the down side, the realization of k — lo diagram having 
a dispersion curve with DBE, requires more sophisticated periodic arrays, such as the 
one shown in Fig. [SJ 

The rest of the paper is organized as follows. 

In Section 2 we briefly outline the electrodynamics of periodic stratified media 
and examine the relation between anisotropy of the layers and the k — oj diagram of 
the stack. We show that at the frequency of degenerate band edge d, the 4x4 transfer 
matrix of a unit cell L cannot be diagonalized or even reduced to the block-diagonal 
form. Based on this, we establish necessary symmetry conditions for a periodic 
stack to develop a degenerate band edge JBJ) and to display the peculiar resonance 
properties associated with it. We prove that such periodic stacks must have at least 
two misaligned anisotropic layers in a unit cell, as shown in the example in Fig. |3J 

In Section 3, we consider the scattering problem for a finite periodic stack. 
We analyze the eigenmode composition of electromagnetic field at the frequency of 
transmission resonance near degenerate band edge. We show that in contrast to the 
case (0 of a regular band edge, in the degenerate band edge case (J7J) the resonance 
field inside the stack does not reduce to a superposition of forward and backward 
propagating cigenmodes. Instead, the contribution of evanescent eigenmodes becomes 
equally important and leads to much stronger dependence (JZJ) of the resonance field 
intensity on the number of layers in the stack. 

The physical and geometrical parameters of stacks used for numerical simulations 
are specified in the Appendix. 
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2. Electrodynamics of periodic stacks of anisotropic layers 

This section starts with a brief description of some basic electrodynamic properties of 
periodic layered media composed of lossless anisotropic layers. Then we turn to the 
particular case of periodic stacks with degenerate band edge. The scattering problem 
for periodic finite stacks, including the Fabry-Perot cavity resonance in the vicinity of 
degenerate photonic band edge will be considered in the next section. 



2.1. Transverse electromagnetic waves in stratified media 

Our consideration is based on time-harmonic Maxwell equations in heterogeneous 
nonconducting media 



V x E (r) = i-B (r) , V x H (r) 

c 



-i-D(r). 



(8) 



where electric and magnetic fields and inductions are related by linear constitutive 
equations 

D (r) = g [f) E (r) , B (f) = £(r)H (f) . (9) 
In further consideration we assume that: 

(i) The direction of plane wave propagation coincide with the normal z to the layers. 

(ii) The second rank tensors e (f) and fi (r) are dependent on a single Cartesian 
coordinate z, normal to the layers. 

(iii) The z direction is a two-fold symmetry axis of the stack, implying that the 
anisotropy axes of individual layers are either parallel, or perpendicular the z 
direction. The latter defines the case of in-plane anisotropy. 

Under the above restrictions, the normal field components E z and H z of 
electromagnetic wave are zeroes, and the system © of six time-harmonic Maxwell 
equations reduces to the following system of four ordinary linear differential equations 
for the transverse field components 



-^-^ (z) = i — M (z) ^ (z) , where * (£) = 

oz c 



E x (z) 

Ey(z) 

H x {z) 

Hy(z) 



(10) 



The 4x4 matrix M (z) in l|10|) is referred to as the (reduced) Maxwell operator. 

In a lossless nonmagnetic medium with in-plane anisotropy, the electric 
permittivity and magnetic permeability tensors have the form 



£ = 







-xy 

s yy 




, A = i- (ii) 

This yields the following explicit expression for the Maxwell operator M (z) in (|10() 








M(z) 












1 " 








-1 







~ £ yy 



















(12) 



where the components of the permittivity tensor may vary from layer to layer. 



Gigantic transmission band edge resonance in periodic stacks of anisotropic layers 10 



2.1.1. The transfer matrix formalism The Cauchy problem 

for the reduced Maxwell equation QlOjl has a unique solution 

*(z)=T(z,zo)*(zo), 



=i~M(z)*(z), 9[zo) 
oz c 



(13) 



(14) 



where the 4x4 matrix T (z, zq) is referred to as the transfer matrix. The transfer 
matrix l|14|) uniquely relates the values of time-harmonic electromagnetic field * at 
any two points z and zq of the stratified medium. From the definition (|14f) . it follows 
that 

T (z, zo) = T (z, z') T (z', zo) , T (z, z ) = T~ x (z , z) , T (z, z) = I. (15) 
The transfer matrix of a stack of layers is defined as 

T S = T(D,0), 

where z = and z — D are the stack boundaries. The greatest advantage of the 
transfer matrix formalism stems from the fact that the transfer matrix of an arbitrary 
stack is a sequential product of the transfer matrices T m of the constitutive layers 



T s = l[T n 



(16) 



If the individual layers m are homogeneous, the corresponding single-layer transfer 
matrices T m can be explicitly expressed in terms of the respective Maxwell operators 
M m 

T m = exp (iD m M m ) , (17) 

where D m is the thickness of the m-th layer. The explicit expression for the Maxwell 
operator M m of a uniform dielectric layer with in-plane anisotropy is given by Eq. i|12|) . 
Thus, Eq. 116(1 together with (|17|) and (|12f> provide an explicit analytical expression 
for the transfer matrix T$ of a stack of dielectric layers with in-plane anisotropy. 

The 4x4 transfer matrix of an arbitrary lossless stratified medium displays the 
fundamental property of J-unitarity ^H] 



T f = JT- 1 J, where J 





0-1 
0-10 

1 



(18) 



which implies, in particular, that 



IdetTl = 1. 



(19) 



Different versions of the 4x4 transfer matrix formalism have been used in 
electrodynamics of stratified media composed of birefringent and/or gyrotropic layers 
for decades (see, for example, ^3 El ID] an d references therein). In this paper 
we use exactly the same notations and terminology as in our previous publications 
[131 1141 1151 IT^] on electrodynamics of stratified media. 
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2.2. Eigenmodes in periodic layered media 

In a periodic layered medium, all material tensors are periodic functions of z, and so 
is the 4x4 matrix M(z) in IjlUfl . Usually, the solutions ^ k (z) of the reduced Maxwell 
equation (f 1 01) with the periodic M(z) can be chosen in the Bloch form 

* fe (z + L) = e lkL V k (z) , (20) 

where the Bloch wave number k is defined up to a multiple of 2ir/L. The definition 
Ijl4|l of the transfer matrix together with Eq. 1)20(1 yield 

T(z + L lZ )^ k (z) = e lkL ^ k (z). (21) 

Introducing the transfer matrix of a unit cell L 

T L = T(L,0), (22) 

we have from Eq. I|21[) 

T L $ fe - e 4fei $ fc , where $ fc = * fc (0) . (23) 

Thus, the eigenvectors of the transfer matrix Tl of the unit cell are uniquely related 
to the Bloch solutions \&fc (z) of the reduced Maxwell equation (|10|l 

$ 4 = ^(0), i= 1,2,3,4. (24) 

The respective four eigenvalues 

X, = e lk ' L , i = 1,2,3,4 (25) 

of Tl are the roots of the characteristic polynomial F4 (X) of the forth degree 

F 4 (X) = dot {T L - XI) = 0. (26) 

Unit cell of a periodic stack can be chosen differently. For example, the choice 
A\ — B — A2 specified in Fig. 0is as good as A2 — A\ — B. Different choice of a unit cell 
corresponds to its shift in the z direction and results in the following transformation 
of the respective transfer matrix Tl 

T' L =T (0, Z) T L T (Z, 0) = [T (Z, O)]" 1 T L T (Z, 0) , (27) 

where Z is the amount of the shift. The modified transfer matrix is similar to the 
original one and has the same set of eigenvalues. 

For any given u>, the characteristic equation (|26(l defines a set of four eigenvalues 
J2HJ)- Real k (or, equivalently, \X\ — 1) correspond to propagating Bloch modes, while 
complex k (or, equivalently, \X\ ^ 1) correspond to evanescent modes. 

The J-unitarity i|18|) of Tl imposes the following restriction on its eigenvalues 

{Xr 1 } ={X*}, * = 1,2,3,4 (28) 

or, equivalently 

{h} = {k*}, i = 1,2,3,4, (29) 

for any given uj. In view of the relations (|29|l or (|28|l . one can distinguish the following 
three different situations. 
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A. All four wave numbers are real 

k\ = ki, k,2 = k$ = fcg, k/± = k\. (30) 
In the example in Fig. H3 this case relates to the frequency range 

0<Lu<Lu a . (31) 

In this case, all four Bloch eigenmodes are propagating. 

B. Two wave numbers are real and the other two are complex 

k\—k\, k2 — k%, ki = fcj, where ks ^ k%, k/±^k\. (32) 
This case relates to the frequency range 

LU a < LU < LU g , (33) 

in Fig. Eta), or the frequency range 

LU a < LU < LU dl (34) 

in Fig. HJb). In both cases ll^.'iH and (|34|l . two of the four Bloch eigenmodes are 
propagating and the remaining two are evanescent with complex conjugated wave 
numbers. 

C. All four wave numbers are complex 

&2 = k*i &4 = fcg, where k± 5^ k*, ^ fe^ / ^3, k^^k\. (35) 

This situation relates to a frequency gap, where all four Bloch eigenmodes are 
evanescent. In the example in Fig. El this case corresponds to 

LOg < LU, Or LUd < LU. (36) 

Notice that the case (B) of two propagating and two evanescent modes can only 
occur in periodic stacks of anisotropic layers. If all the layers in a unit cell are isotropic, 
the four Bloch eigenmodes are either all propagating or all evanescent, as is the case 
in Fig. ^ For a given frequency lu, the four Bloch eigenmodes correspond to two 
different polarizations and two opposite directions of propagation. 

2.2.1. Non-Bloch solutions at stationary points of the k — lu diagram So far, we 
have considered only the cases where all four solutions for the reduced Maxwell 
equation (|10f> can be chosen in the Bloch form H20|) . All such cases fall into one 
of the following three categories: (A) all four eigenmodes are propagating, (B) two 
modes are propagating and the other two are evanescent, (C) all four eigenmodes are 
evanescent. This classification of the eigenmodes does not apply at the frequencies of 
stationary points on the k — lu diagram, where the group velocity u of some of the 
propagating modes vanishes 

u = dLu/dk = 0. (37) 

At a stationary point of a dispersion curve, not all four solutions of the Maxwell 
equation Hl()[) are Bloch waves, as defined in H2()(l . Instead, some of the solutions can 
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EJ Using these 
HJa), the four 



be algebraically diverging non-Bloch eigenmodes. Such eigenmodes can be essential 
for understanding the resonance effects in finite and semi-infinite periodic arrays. 

Let us consider the situation at stationary points l|37fl on the k — u> diagram 
in terms of the matrix Tl- Although at any given frequency oj, the reduced Maxwell 
equation lfTU|) has exactly four linearly independent solutions, it does not imply that the 
respective transfer matrix Tl in (|23[1 must have four linearly independent eigenvectors 
J23J. Indeed, although the matrix Tl is invertible, it is neither Hermitian, nor unitary 
and, therefore, may not be diagonalizable. Specifically, if the frequency approaches 
one of the stationary points (|37[) . some of the eigenvectors in (|23[) become nearly 
parallel to each other. Eventually, as u> reaches the the stationary point value, the 
number of linearly independent eigenvectors <5>k becomes lesser than four, and the 
relation (|24J) does not apply at that point. 

Examples of different stationary points 1|37[) are shown in Fig. 
examples, let us take a closer look at these special frequencies. 

At the frequencies oj g of the photonic band edge g in Fig. 
eigenmodes include: 

- one propagating mode with k = ir/L and zero group velocity, 

- one non-Bloch eigenmode linearly diverging with z, 

- a pair of evanescent modes with equal and opposite imaginary wave numbers. 

At the frequencies uj a corresponding to the point a in Fig. El the four solutions 
of Eq. (jTUjl include: 

- the propagating mode with k = tt/L and zero group velocity, 

- one non-Bloch eigenmode linearly diverging with z, 

- a pair of propagating modes having equal and opposite group velocities and 
belonging to the dispersion curve other than the one containing the point a. 

Of special interest here is the frequency u>d of degenerate photonic band edge d 
in Fig. EJb). In this case, the four solutions of Eq. Ijl0|l include 17 : 

- the propagating mode with fc = n/L and zero group velocity, 

- three non-Bloch eigenmodes diverging as z, z 2 , and z 3 , respectively. 

At any particular frequency, the existence of non-Bloch eigenmodes can be directly 
linked to the canonical Jordan form of the respective transfer matrix Tl ■ Indeed, the 
4x4 matrix Tl, being invertible, can have one of the following five different canonical 
forms 



Ti = 



T 2 i = 



T, = 
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of which all but T\ have non-trivial Jordan blocks and, therefore, are not 
diagonalizable. Each m x m Jordan block is associated with a single eigenvector 
of the respective T - matrix. Therefore, the total number of eigenvectors of different 
4x4 matrix in (|38[1 is lesser than the number four of the solutions for the reduced 
Maxwell equation (|10|l . The only exception is the diagonalizable case X\, where the 
four Tl eigenvectors correspond to four Bloch eigenmodes, as prescribed by l|24(l . 
Generally, each nontrivial m x m Jordan block of the T - matrix is associated with 
the m eigenmodes of Eq. I)10[l. of which one is a propagating Bloch eigenmode with 
zero group velocity, and the other m — 1 are non-Bloch eigenmodes algebraically 
diverging with z. The details can be found in any course on linear algebra, for example, 
^OEOEnil^- Let us cons id er each case separately. 

At general frequencies, different from those of stationary points \XJ\. the 4x4 
matrix Tl is always diagonalizable. Its canonical (diagonalized) form is trivial and 
coincides with T\ in i|38|) . The four eigenvectors are defined in (|24H . All the possibilities 
here reduce to one of the three cases (A), (B), or (C) described earlier in this Section. 
None of them involves non-Bloch solutions. 

At the frequency uj g of a regular photonic band edge in Fig. El^a), the canonical 
Jordan form of the respective transfer matrix is T21 in (|38|l . where 

x x = -1, x 3 = x; = x- 1 ? 1. 

The 2x2 Jordan block relates to one propagating Bloch mode with k = ir/L and 
zero group velocity, and one non-Bloch linearly diverging eigenmode. The pair of real 
eigenvalues A3 and X4 — X^ 1 relate to the pair of evanescent modes at uu = uo g . 

By contrast, at the frequency of photonic band edge g in Fig. ^b) related to the 
stack of isotropic layers, there is no evanescent modes. The canonical Jordan form of 
the respective transfer matrix coincides with T22 in {HHt ■ where 

Xi = X 2 = -1. 

Each of the two 2x2 identical Jordan blocks relates to one propagating Bloch mode 
with k = n/L and zero group velocity, as well as one non-Bloch linearly diverging 
eigenmode. The two Jordan blocks of T22 correspond to two different polarizations of 
light. 

At the frequency oj a in Fig. the canonical Jordan form of the respective transfer 
matrix is T21, where 

X 1 = -l, X 3 =Xl |A 3 | = |A 4 | = 1. 

The double eigenvalue X\ — — 1 of the 2x2 Jordan block relates to one propagating 
Bloch mode with k — n/L and zero group velocity, and one non-Bloch linearly 
diverging eigenmode. The pair of complex eigenvalues X 3 and X4 = X 3 relate to the 
pair of propagating modes having equal and opposite group velocities and belonging 
to the dispersion curve other than the one containing the point a. 

The canonical form T3 in (|38f) relates to a k — uj diagram with stationary inflection 
point. Such a stationary point cannot be realized in a reciprocal periodic stack at 
normal light propagation |14l lift) . 

Of particular interest here is the case JHJ of degenerate band edge. At the 
degenerate band edge frequency ojd, the canonical Jordan form of the respective 
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transfer matrix is T 4 in l|38|) 



Tr ~ Ta 



X 1 

X 









1 

X 



X 



where X = ±1 



(39) 



More specifically, in the case shown in Fig. HJb) 

X = X rl = e ikdL = -1, at u = Lu d . 



d = z a 

The matrix l|39l) presents a single 4x4 Jordan block and has a single eigenvector, 
corresponding to the propagating eigenmode with k = n/L and zero group velocity. 
The other three solutions for the Maxwell equation (|10|) at u> — ojd are non-Bloch 
eigenmodes diverging as z, z 2 , and z 3 , respectively. 

If the frequency to deviates from the stationary point (|37|l . the transfer matrix Tj, 
become diagonalizable with the canonical Jordan form T\ in H38J) . The perturbation 
theory relating the non-Bloch eigenmodes at the frequency of degenerate band edge 
to the Bloch eigenmodes in the vicinity of this point is presented in |17j . 

2.3. Symmetry conditions for the existence of degenerate band edge 

Not any periodic stack can develop a degenerate band edge, defined in ijfJJl. Some 
fundamental restrictions can be derived from symmetry considerations. These 
restrictions stem from the fact that at the frequency ujd of degenerate band edge, 
the transfer matrix Tl must have the Jordan canonical form (|39[1 . Such a matrix 
cannot be reduced to a block-diagonal form, let alone diagonalized. Therefore, 

- a necessary condition for the existence of degenerate band edge is that the 
symmetry of the periodic array does not impose the reducibility of the transfer matrix 
Tl to a block- diagonal form. 

The above condition does not imply that the transfer matrix Tl must not be 
reducible to a block-diagonal form at any frequency uj on the k — uj diagram. Indeed, 
at a general frequency w, the matrix Tl is certainly reducible, and even diagonalizable. 
The strength of the symmetry imposed reducibility is that it leaves no room for 
exceptions, such as the frequency ujd of degenerate band edge, where the transfer 
matrix Tl must not be reducible to a block-diagonal form. Therefore, in the case of 
symmetry imposed reducibility, the very existence of degenerate band edge is ruled out. 
Observe that in most periodic layered structures, the symmetry of the periodic array 
does require the matrix Tl to be similar to a block-diagonal matrix at all frequencies. 
In all these cases, the stack symmetry is incompatible with the existence of degenerate 
band edge on the k — uj diagram. 

Let us apply the above criterion to some specific cases. 

In periodic stacks of isotropic layers, the Maxwell equations for the waves with the 
x- and the y - polarizations are identical and decoupled, implying that the respective 
transfer matrix can be reduced to the block-diagonal form 



Tr 





Tl2 








T 2i 


T22 














T11 


T12 








T21 


T22 



(40) 
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The two identical blocks in (|4L)fl correspond to two different polarizations of light. 
The characteristic polynomial F^X) of the block-diagonal matrix l|4(J|) factorizes into 
a product of two identical second degree polynomials related to electromagnetic waves 
with the x- and the y - polarizations, respectively 



Fi(X) = F 2 {X)F 2 (X). 



(41) 



The block-diagonal structure of the matrix (|40|l rules out the existence of degenerate 
band edge in periodic stacks of isotropic layers. In fact, the transfer matrix Tj, in this 
case can only have the following two canonical forms 



Ti 



X 














x- 1 














X 














x- 1 



To. 



±1 


1 











±1 














±1 


1 











±1 



compatible with (|40|l . The case T\ relates to a general frequency, while the case T22 
relates to a photonic band edge, like the one shown in Fig. ^b). 

Let us now turn to the situation where all or some of the layers of the periodic 
stack are birefringent. The in-plane dielectric anisotropy (|11J) may allow for degenerate 
band edge on the k — uj diagram, but not automatically. 

Let us start with the simplest periodic array in which all anisotropic layers of the 
stack have aligned in-plane anisotropy. The term " aligned" means that one can choose 
the directions of the in-plane Cartesian axes x and y so that the permittivity tensors 
in all layers are diagonalized simultaneously. In this setting, the Maxwell equations 
for the waves with the x- and the y - polarizations are still separated, implying that 
the respective transfer matrix can be reduced to the block-diagonal form 



Tr = 



Tn 


Tl 2 








T21 


T22 














T33 


T34 








T43 


T44 



(42) 



The two blocks in l|42l) correspond to the x and y polarization of light. The forth 
degree characteristic polynomial of the block-diagonal matrix (|42|l factorizes into the 
product 

F 4 (X) =F x (X)F y {X), (43) 

where F X (X) and F y (X) are independent second degree polynomials related to 
electromagnetic waves with the x- and the y - polarizations, respectively. The typical 
k — lo diagram in this case will be similar to that shown in Fig. I^a) with two separate 
curves related to two linear polarizations of light. Again, the block-diagonal structure 
of the matrix Q42[) rules out the existence of degenerate band edge in periodic stacks 
with aligned anisotropic layers. The transfer matrix Tl in this case can only have the 
following two canonical forms 



Ti = 



x x 




























X 2 














at 1 



; T 2 2 — 



±110 
±1 
0X0 

000 x- 1 
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compatible with |@2J|. The case T\ relates to a general frequency, while the case T22 
relates to a photonic band edge. 

As we have seen, the presence of anisotropic layers may not necessarily lift 
the symmetry prohibition for the degenerate band edge, because the symmetry of 
the periodic array may still be incompatible with the canonical Jordan form l|39|) . 
Generally if the space symmetry group G of the layered structure includes a mirror 
plane m| parallel to the z direction, this would guarantee the reducibility of the matrix 
Tl to a block-diagonal form. Indeed, a standard line of reasoning gives that if mu G G, 
and the y axis is chosen perpendicular to the mirror plane mu, then the waves with 
the x- and the y - polarizations have different parity with respect to the symmetry 
operation of reflection and, therefore, are decoupled. The latter leads to reducibility 
of the matrix Tl to the block-diagonal form (|42(l . Thus, a formal necessary condition 
for the existence of degenerate band edge on the k — lu diagram can be written as 
follows 

m|| i G. (44) 

None of the common periodic layered structures satisfies this criterion and, 
therefore, none of them can develop the degenerate band edge. For example, even 
if anisotropic layers are present, but the anisotropy axes in all anisotropic layers are 
either aligned, or perpendicular to each other, the symmetry group G of the stack 
still has the mirror plane m||, which guarantees the separation of the x- and the y - 
polarizations and the reducibility of the respective transfer matrix Tl to the block- 
diagonal form H42J1 . The only way to satisfy the condition (|44J) and, thereby, to allow 
for degenerate band edge on the k — lu diagram, is to have at least two misaligned 
anisotropic layers in a unit cell with the misalignment angle being different from 
and 7r/2, as shown in the example in Fig. [S] 

Observe that the presence of B layers in the periodic array in Fig. [S] is also 
essential, unless the two anisotropic layers A\ and A2 have different thicknesses or are 
made of different anisotropic materials. In Fig. the layers Ai and A2 differ only 
by their orientation in the x — y plane, but otherwise, they are identical. In such a 
case, if the B layers are removed, the point symmetry group of the periodic stack in 
Fig. [3] rises from D2 to -D2/1 acquiring the glide mirror plane mu. This, according 
to the criterion (|44() . imposes the reducibility of the matrix Tl to the block-diagonal 
form H42JI . regardless of the misalignment angle between the adjacent A layers. The 
symmetry imposed reducibility rules out the possibility of the degenerate band edge 

The k — lu diagram of the periodic stack in Fig. [S]with the B layers removed is 
shown in Fig. E>[d). 

In the numerical example considered in the next section, the B layers are simply 
empty gaps of certain thickness Db between the adjacent double layers A\— A2- The 
misalignment angle is chosen ir/4. By changing the thickness Db of the gap B, one 
can change the k — lu diagram of the periodic stack, as shown in Fig. El Similar effect 
can be achieved by changing the misalignment angle between the adjacent A layers. 

3. Transmission resonance in the vicinity of degenerate band edge 

3.1. Scattering problem for periodic semi-infinite stack 

To solve the scattering problem for a plane monochromatic wave incident on a finite 
stack of anisotropic layers we use the following standard approach based on the 4x4 
transfer matrix. 
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Let (z), \I/p (z), and (z), be the incident, reflected, and passed plane waves 
in vacuum. Allowing for general elliptic polarization of the waves, we have 



Rx 



R 



*p(*) = 



P, 
Py 
~Py 

Pa-. 



The respective domains of definition are 



z < 0, for (z) and *p (z) , 
D < z, for *p (z) . 



(45) 



where D is the stack thickness. The field inside the stack is denoted by *S?t {z). Except 
for the stationary points (|37|l on the k — u> diagram, (z) can be decomposed into 
a superposition of the four Bloch solutions i|2U[l of the Maxwell equation I|1U|) 



(46) 



This representation is meaningful if the periodic stack contains a significant number 
of unit cells L. Otherwise, if there are just a few layers in the stack, the representation 
H46fl is formally valid, but not particularly useful. 

The boundary conditions at the two slab/ vacuum interfaces are 



* (0) = tf , (0) + * R (0) ,*(D) = * P (D) 
The transfer matrix Ts of a periodic stack is 



T s = (T L ) 



(47) 



(48) 



The relation 



* (£>) = (0) , (49) 

together with the pair of boundary conditions l|47|) allow to express both the 
reflected wave "Jp and the wave ^p passed through the slab, in terms of a given 
incident wave ^1 and the elements of the transfer matrix Ts- This also gives the 
transmittance/reflectance coefficients of the slab defined as 



td 



Sp 
Si 



\*p(D)\ 
I*/ (0)1 



2~) PD — ~ 



Si |*/(0)| 2 



(50) 



where S = cW is the Poynting vector of the respective wave. In the case of a lossless 
stack 

td +P D = !• 

The field distribution (z) inside the slab is found using either of the following 
expressions 

9 T {z)=T(z,0)[*i(0) + * a (0)]=T(z,D)* P (D), < z < D. (51) 



The above procedure is commonly used for the frequency-domain analysis of 
periodic and non-periodic layered structures involving anisotropic and/or gyrotropic 
layers (see, for example, QUI i i EH 1 > EH El E3 EH an d references therein). 
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In addition to the field distribution inside the slab, we are also interested in 
its eigenmode composition. The latter is particularly important since it allows to 
explain the fundamental difference between Fabry-Perot resonance in the vicinity of a 
degenerate band edge and a similar resonance in the vicinity of a regular band edge. 
Throughout this section we consider only the first transmission resonance, closest to 
the respective band edge. 

3.2. Field composition at the frequency of transmission resonance 

In the case of transmission resonance in a periodic stack of isotropic layers, the 
resonance held inside the stack is a simple standing wave composed of two propagating 
Bloch modes with opposite group velocity (see Eqs. through (J5J and the comments 
therein). The introduction of anisotropy in itself does not change qualitatively 
the resonance picture, as shown in Fig. The only difference is that the 
electromagnetic field "J/t (z) inside the stack can now have both propagating and 
evanescent components. But at the frequency of a transmission resonance, the 
contribution of the evanescent components is negligible, as shown in Fig. 1101 So, 
basically, one can still see the resonance field inside the stack as a simple standing 
wave composed of a pair of propagating modes with greatly enhanced amplitude, 
compared to that of the incident wave. The formulas (JTJ through (JSJ still apply here, 
provided that the number N of unit cells in the stack is not too small. 




Figure 9. Transmission resonance in the vicinity of regular photonic band edge 
(BE) g in Fig. 6(a) in the stack composed of 32 unit cells L. a) Transmission 
dispersion at the frequency range including the regular BE g and two closest 
transmission resonances, b) Smoothed field intensity distribution (59) inside the 
stack at the frequency ui\ of the first transmission resonance 1. 



In Fig. HUf a) we show the squared amplitudes (intensity distributions) 



\*j(z)\ 2 , j = 1,2,3,4, (52) 



of the individual Bloch components of the resulting field *S?t (z) at the frequency of the 
first transmission resonance near the regular band edge g in Fig. E{a). The numbers 
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1 and 2 designate the forward and backward propagating components of the standing 
wave. The evanescent contributions 3 and 4 are negligible. In Fig. HOf b) we show 
the squared amplitudes of the combined contribution of the pair of propagating waves 
(pr), and the combined contribution of the pair of evanescent waves (ev) 

IV (z)| 2 = |*i (z) + * 2 (z)| 2 , \* ev (z)\ 2 = 1*3 (z) + * 4 (z)| 2 . (53) 

Obviously, 

* T (z) « (z) . (54) 




Figure 10. The Bloch composition of the resonant field in Fig. 9(b) 
corresponding to the transmission peak 1 in Fig. 9(a). a) Square moduli of each 
of the four Bloch components; 1 and 2 are forward and backward propagating 
components, b) Square moduli of the combined contribution of two propagating 
components (pr), and the combined contribution of two evanescent components 
(ev). The evanescent modes contribution is negligible. 

Now let us turn to the case of transmission resonance in the vicinity of degenerate 
band edge d in Fig. EJb). There are several features that sharply distinguish this case 
from the similar transmission resonance near the regular band edge g in Fig. |B[a). 

First of all, for a given number N of unit cells in the stack (N = 32), the field 
intensity in Fig. 02b) is by several orders of magnitude higher, compared to that in 
Fig. Etb), in spite of the close similarity of all numerical parameters of the respective 
periodic structures (see the Appendix). In the case of transmission resonance in the 
vicinity of degenerate band edge (DBE), the resonance field intensity increases as TV 4 , 
while in the case of a regular band edge (BE), the field intensity is proportional to 
iV 2 . 

The second distinction is that in the case of degenerate band edge (DBE), the 
field intensity near the slab boundaries at z = and z = D, increases as 

DBE case: |* T (z)| 2 oc z 4 , (D — zf . (55) 

By contrast, in the case of a regular band edge (BE), the field intensity near the stack 
boundaries rises at a much slower rate 



Regular BE case: |* T (z)| 2 oc z 2 , (D - zf , 



(56) 
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which is characteristic of a regular standing wave composed of two propagating 
components. 

Finally, in the case of a regular band edge, the transmission resonance field i z ) 
is a standing wave composed of two propagating Bloch modes with opposite group 
velocities. The evanescent modes do not participate in the formation of the resonance 
field, as clearly seen in Fig. By contrast, in the case of transmission resonance 
in the vicinity of DBE, the role of evanescent components in the formation of the 
resonance field is absolutely crucial. Indeed, as shown in Fig. El amplitude of 
the propagating and evanescent components are comparable in magnitude. More 
importantly, the combined contribution ^> pr (z) of the two propagating components 
does not even resemble a standing wave with the nodes at the slab boundary, as was 
the case in Fig. fTHT b). Instead, comparing Fig. E{b) and Fig. fTTT bl we see that at 
the slab boundaries at z = and z = D. the propagating and evanescent components 
interfere destructively, almost canceling each other 

at z = and z = D: ^ pr (z) w —^ e v (z) , 

while the individual Bloch components (|52l) remain huge. We remind that in all cases, 
the intensity of the incident wave ^/ is unity. 




Figure 11. The Bloch composition of the resonant field in Fig. 8(b) 
corresponding to the transmission peak 1 in Fig. 7. a) Square moduli of each 
of the four Bloch components, 1 and 2 are forward and backward propagating 
components, b) Square moduli of the combined contribution of two propagating 
components (pr), and the combined contribution of two evanescent components 
(ev). Propagating and evanescent contributions interfere destructively at the slab 
boundaries. 

Qualitatively, such a bizarre resonance behavior in the vicinity of the degenerate 
band edge can be characterized as follows. As we already mentioned, at the frequency 
uod of degenerate band edge, three of the four solutions for the Maxwell equation 
<|1(J|I are non-Bloch eigenmodes diverging as z, z 2 , z 3 . Although the frequency w\ 
of the transmission cavity resonance is slightly different from and, therefore, all 
four eigenmodes are Bloch waves, the close proximity of lj\ to u>d causes all the 
abnormalities seen in Fig. 1101 To describe these behavior in mathematically consistent 
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way, one should start with the one solutions at u> = to,} as zero approximation. Then 
using the perturbation theory for the non-diagonalizable transfer matrix Tj, (u>d), one 
can derive an asymptotic theory for the case of large N. The perturbation theory 
for the degenerate band edge d was developed in ^7]. It turns out that near the slab 
boundaries at z -C D and (D — z) <§C D, the field ^>t (z) at the transmission resonance 
frequency w\ is well approximated by quadratically diverging non-Bloch eigenmode, 
corresponding to u — u>d- It qualitatively explains an extremely rapid growth l|55|l of 
the resonance field inside the slab as one moves away from either slab boundary. It 
also explains the sharp dependence J7J) of the resonance field intensity on the number 
N of unit cells in the stack. 
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4. Appendix 



Numerical parameters of layered arrays 

The periodic array in Fig. [S] has three layers in a unit cell L, of which two {A\ and 
A2) are anisotropic and have the same thickness Da- The third layer B is isotropic 
with the thickness Db- In our numerical simulations, the B layers are empty gaps of 
variable thickness. 

The anisotropic layers A 1 and A 2 are made of the same lossless dielectric material 
and have the same thickness. The respective dielectric permittivity tensor is 

Sa + S cos 2ip S sin 2<p 
S sin 2tp sa — 5 cos 2(p 











where S describes the magnitude of in-plane anisotropy, while the angle ip defines the 
orientation of the anisotropy axes of the respective layer in the x — y plane. The most 
critical parameter of the periodic structure in Fig. [3] is the misalignment angle 

f = <Pi - ¥2 (57) 

between the adjacent anisotropic layers A\ and A 2 . This angle determines the 
symmetry of the periodic array and, eventually, what kind of k — w diagram it can 
display. It is important for our purposes that the misalignment angle is different from 
and 7r/2. In our numerical simulations we set 

<Pi = 0)^2 = V 4 , 

and use the following expressions for the dielectric permittivity of the A layers 
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In our numerical simulations we set 
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(58) 



sa = 13. 61, S 

The numerical value of e zz is irrelevant. 

In the case of the periodic array of isotropic layers in Fig. (Ha), we set 

e A = 3.78, 6 = 0, D A = D B = 0.5 x L. 



4-2. Description of plots 

In all plots of the field intensity distribution we, in fact, plotted the following physical 
quantity 

(|* (z)\ 2 ) = (E (z) ■ E* (z) + H(z)- H* (z)) l , (59) 

which is the squared field amplitude averaged over local unit cell. The real 
electromagnetic energy density W (z) is similar to |$ (z)\ 2 . Both of them are strongly 
oscillating functions of the coordinate z with the period of oscillations coinciding with 
the unit cell length L. Thus, the quantity (|59"|l can be interpreted as the smoothed 
field intensity distribution, with the correction coefficient of the order of unity. 

In all plots, the wave number k and the frequency u> are expressed in units of 
and cL _1 , respectively. 



